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PREFACE 

This  Memorandum  la  part  of  RAND's  study  for  the  Advanced  Research 
Projects  Agency  of  those  phenomena  which  affect  the  performance  of 
optical  reconnaissance  and  guidance  equipment.  The  objective  of  these 
studies  is  to  provide  sufficient  understanding  to  permit  the  systems 
analyst  to  compute  performance  estimates  under  various  operational 
conditions. 

In  optical  communications  and  related  devices,  the  random  vari¬ 
ations  in  the  received  signal  due  to  atmospheric  turbulence  can 
represent  a  severe  limitation  to  system  performance.  This  Memorandum 
gives  an  improved  expression  for  the  effect  of  a  finite  diameter 
receiving  aperture  in  reducing  the  variance  of  a  fluctuating  light 
signal,  and  should  be  of  interest  to  those  concerned  with  the  use  of 


lasers  in  the  atmosphere. 
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SUMMARY 

The  aperture -averaging  faetar  of  a  circular  aparture  is  given. 

This  factor  givaa  tha  effect  of  a  finite  receiving  aperture  on  spherical 
and  plane  waves  in  reducing  the  variance  of  a  fluctuating  light  signal. 
Curves  of  the  reduction  factor  and  normalized  signal  standard  deviation 
ae  a  function  of  range  and  receiver  aperture  diameter  are  presented  and 
are  compared  with  thoee  which  were  previously  calculated  by  Fried. 

It  is  shiwn  that  Fried’s  expression  agrees  with  the  results  obtained 
here  only  for  propagation  distances  much  less  than  the  propagation 

distance  where  the  average  field  is  down  by  a  factor  of  the  order 

-1 
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I.  INTRODUCTION 

When  a  large  receiving  aperture  la  ueed  to  collect  scintillating 
light,  such  as  laser  light  that  has  propagated  through  a  turbulent 
medium,  the  percentage  fluctuation  is  not  as  large  as  would  be  obtained 
If  a  small  aperture  were  used.  The  aperture  performs  an  averaging 
process.  This  Is  explained  by  the  fact  that  for  D  >  pQ,  where  pQ  Is 
the  characteristic  length  of  the  log-amplitude  covariance,  several 
field  inhomugfcui.x^i.t.j  with  uitierenL  signs  can  oe  lound  within  the 
limits  of  the  receiving  aperture,  and  therefore  they  partially  compensate 
one  another.  In  the  inertial  subrange,  pQ  =- /Tz,  where  D  is  the  diameter 
of  the  assumed  circular  receiving  aperture,  \  is  the  optical  wavelength, 
and  z  is  the  propagation  distance. 

This  process  has  been  studied  by  Fried, ^  who  calculates  the 

aperture -averaging  factor  (defined  below)  for  a  large  range  of  the 

relevant  parameters.  Central  to  his  calculation  is  the  covariance  of 

irradiance,  which  he  takes  as  that  result  that  would  be  obtained  from 

the  Rytov  approximation.  As  has  been  discussed  elsewhere, the 

Rytov  approximation  is  invalid  for  propagation  distance  z  exceeding  a 

critical  distance  z^.  The  range  z^  is  the  range  at  which  the  average 

component  of  the  field  is  down  by  a  factor  of  the  order  e  ^  and  is 

typically  ~  1  km  for  horizontal  propagation  paths  near  the  earth's 

surface.  The  Rytov  approximation  predicts  large  dispersion  of  the 

intensity  for  z  >  zc>  whereas  experimentally/ 3^  the  dispersion  of 

the  intensity  tends  asymptotically  to  unity  for  z  >  z  .  Hence 

calculations  based  on  the  Rytov  approximation  are  invalid  for  z  >  z  . 

c 

Here  we  derive  the  aperture -averaging  factor  based  upon  a  functional 
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form  of  the  dispersion  of  the  intensity  that  is  consistent  with  the 
experimental  observations.  Both  plane  and  spherical  waves  are  dis¬ 
cussed.  The  effect  of  a  finite  transmitting  aperture  is  not  considered 
here.  Only  propagation  paths  with  constant  turbulence  parameters  are 
discussed  (e.g.,  horizontal  paths),  since  the  dependence  of  the  tur¬ 
bulence  parameters  on  height  is  not  known  well  and  unduly  complicates 
the  problem  mathematically. 


3 


II.  Amw-AVPACUW  FACTOR 

Consider  a  circular  receiving  aperture  with  diameter  D.  The  signal 

2  (1) 

variance  ag  is  given  by 

a]  -  2it  dpp  K(p  ,D)  Cr(p)  (1) 

where 

(?Vl(§) 

K(p,D)  -  J 

(  0  if  p  >  D 


Cr(p)  -  (l(x)  -  IQ)  (l(x')  -  IQ)  (3) 

IQ  is  the  average  irradiance  in  the  absence  of  turbulence,  p  ■  |x  -  x'|, 
and  the  bar  over  a  quantity  represents  the  ensemble  average  of  the 
quantity.  The  quantity  Cj(p)  is  the  covariance  of  the  Intensity  in 
a  plane  perpendicular  to  the  propagation  direction  at  distance  t  from 
the  source.  The  quantity  K(p,D)  is  a  geometrical  factor  which  is 
seen  to  be  Just  the  area  of  overlap  of  two  circles  of  diameter  D 
whose  centers  are  displaced  a  distance  p. 

2 

For  normalization  purposes,  it  is  convenient  to  divide  a  by 
2  2 

(rrD  /A)  Cj(0)  ,  the  signal  variance  that  would  be  obtained  if  the 

irradiance  fluctuations  were  perfectly  correlated  over  the  entire 
aperture.  We  call  this  ratio  0,  the  aperture-averaging  factor.  We 


-  (§)  [*  -  m 


if  p  £  D 


(2) 


have 
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(A) 


from  which,  by  using  appropriate  values  of  C^Cp),  the  aperture 
averaging  factor  may  be  obtained. 

Fried^^  based  his  calculation  on 


cj(p)  ■  Iq  {exp  CAC^(p)]  -  lj 


(5) 


where  C^(p)  is  the  covariance  of  the  log-amplitude.  Equation  (5)  is 

valid  at  all  ranges  provided  the  log-amplitude  is  normally  distributed. 

Fried  has  computed  the  log-amplitude  covariance  using  a  modified 

Rytov  approximation  (modified  to  be  consistent  with  energy  conser- 

vation)  ,  which  will  be  denoted  by  C^(p).  As  has  been  discussed 

(2) 

elsewhere,  this  approximation  is  invalid  for  z  >  z^.  The  critical 
range  z for  \  ■  0.6328m<  (He  -  Ne  laser  wavelength)  is: 


Structure  Constant 


n2  n  in-15  -2/3 

C  -  0. 3  x  io  cm 
n 


,  v  in-15  -2/3 

3  X  10  cm 


,n  in-15  -2/3 

30  x  10  cm 


Plane  Waves  Spherical  Waves 


z  ■  3. 6  km 

c 


z  ■  5. 9  km 
c 


■1.0km  "1.7km 

■  0.3  km  ■  0.48  km 


No  theoretical  expression  for  C^(p)  which  is  valid  for  arbitrary 
z  exists.  An  expression  which  is  a  smooth  interpolation  between  the 
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limits  of  small  and  large  z  and  gives  a  value  for  the  intensity  vari¬ 
ance  which  agrees  to  within  20  percent  with  the  experimental  values^ 

(2) 

has  been  given  by  Yura: 

CI(p)  “  Io[1  '  exP  ('2  ^(P))]  (6) 

Equation  (6)  agrees  with  Eq.  (5)  for  small  z  (where  cj(0)  «  1)  and 
tends  asymptotically  to  I2  for  large  z  (where  C*(0)  »  1).  Figure 
1  is  a  graph  of  (^(O)/^  versus  C^(0)  as  given  by  Eqs.  (5)  and  (6). 

The  aperture -averaging  factor  0  is  obtained  here  from  (^(p)  as 
given  by  Eq.  (6).  Fried  and  Cloud^  ^  have  given  an  analytic  reprcsen- 
tation  of  C^(p)  for  plane  waves  assuming  horizontal  propagation  paths 
and  that  the  Kolomogorov  spectrum  represents  the  turbulence -induced 
fluctuations  of  the  index  of  refraction.  Their  result  is  given  by 
Eq.  2.12  of  Ref.  4.  Similarly,  Fried^  ^  has  derived  an  analytic 
expression  for  C^(p)  for  spherical  waves.  This  result  is  given  by 
Eq.  2.8  of  Ref.  5. 

We  can  now  calculate  the  normalized  signal  standard  deviation 
l°s/(S)  ]  ,  where  S,  the  average  signal,  equals  ttD2I0/4.  From  Eq.  (4) 
we  find 


The  values  of  C^(0)  have  been  derived  by  Fried  and  Cloud^  for 
plane  waves  and  by  Fried^  for  spherical  waves.  They  are  given  by 
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-L  -7  -2  -1 

10  10  3  10  10  1  10 


c*(0) 

p 

Fig.  1--The  normalized  variance  of  irradiance  C  (0)/K  as  a 

n  I  U 

function  of  log-amplitude  variance  C£(o) 
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C,(0)  ■  0.309k  ^  C2  (plane  waves) 

and 

C*(0>  “  0.124k  z  ^  (spherical  waves) 


(8) 


(9) 


Figure  2  gives  the  aperture -averaging  factor  for  plane  and 

spherical  waves  as  a  function  of  the  normalized  receiver  aperture 

diameter  D // 4z/k  (k  -  2tt/X)  for  z  •  1  km,  X  ■  0.6328m-,  and 
2  -15  -2/3  2 

cn  "  ^  x  0,1  »  where  denotes  the  Index  structure  constant.* 

Figure  3  gives  the  aperture  averaging  factor  for  plane  and  spherical 
waves  as  a  function  of  the  normalized  range  4z/kD2  for  D  -  5  an, 

\  ■  0.6328m-,  and  ■  3  X  10  cm  For  comparison  purposes,  the 

results  based  on  Fried's  expression  are  Included  In  these  graphs. 

The  differences  in  9  for  long  propagation  paths  are  readily  dis¬ 
cernible.  This  is  because  the  aperture  standard  deviation  as  computed 
from  Eq.  (6)  tends  to  a  finite  value  for  large  z,  while  the  same 
quantity  computed  from  Eq.  (5)  increases  without  limit.  Figures  4 
and  5  give  the  aperture-averaging  factor  as  a  function  of  the 
normalized  range  4z/kD2  (for  D  -  5  cm)  and  as  a  function  of  D//4^7k 

(z  -  1  km)  for  X  -  0.6328m-,  and  various  values  of  C2  for  plane  and 

n 

spherical  waves,  respectively.**  Figure  6  is  a  plot  for  plane  waves 


*  2 

The  values  of  Cn  typical  of  turbulence  In  the  first  few  hundred 
meters  of  the  atmosphere  are:  30  x  10-15  m-2/3  (strong  d,yttme  tur_ 
bulence  within  a  few  meters  of  the  ground),  3  x  10‘15  cm"2^  and 
7  x  10-13  cra-Z/3  (moderate  daytime  turbulence  and/or  strong  night¬ 
time  turbulence),  and  0.3  x  10"15  cra"2/3  (very  weak  turbulence  which 
°CCUr#*in  th*  near_neutral  period*  *t  dawn  or  dusk). 

The  specific  expressions  used  here  for  the  normalized  range  and 
receiver  aperture  are  introduced  for  convenience  and  are  not  to  be  In¬ 
terpreted  as  scaling  parameters. 


(x  =  0.6328u,  z  =  1  km,  and  C_  =  3  x  10 


I 
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Fig. 


3-The 


aperture-averaging  factor  e  as  a  function  of 
(x  =  0.6328p,  D  =  5  cm,  and  C2  =  3  x  10"^ 


normalized  range  4z/kD2 


-2/3* 
cm  ) 


D/ywk 


Fig.  4-The  aperture-averaging  factor  e  for  plane  waves  as  a  function 
of  normalized  range  4z/kD^  (A  =  0.6328y,  D  =  5  cm)  and 
normalized  receiver  diameter  D/v^z7k  [\  =  0.6328u, 
z  =  l  km),  for  various  values  of  the  index  structure  constant 
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Fig.  6--The  normalized  signal  standard  deviation  o  for  plane  waves  as  a 
function  of  normalized  range  4z/kD2  (x  =  0.6328u,  D  =  5  cm)  and 
normalized  receiver  diameter  DA/5i?k  (x  =  0.6328y  z  =  1  kml 
for  various  values  of  the  index  structure  constant 
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Of  the  normalized  signal  standard  deviation  c  as  a  function  of  dz/kD2 
and  a,  a  function  of  D/AJ7T  for  the  same  values  of  the  parameters, 
while  Fig.  7  gives  similar  graphs  for  the  normalized  signal  standard 
deviation  for  spherical  waves.  For  comparison  purposes,  the  result, 
of  using  Fried's  expression  are  included  in  Figs.  6  and  7.  At  small 
apertures  the  difference  between  the  present  analysis  and  that  of 
Fried  should  be  readily  observable. 
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Fig.  7--The  normalized  signal  standard  deviation  c  for  spherical  waves 
as  a  function  of  normalized  range  4z/kD?  (x  =  0.6328u,  D  =  5  cm) 
and  normalized  receiver  diameter  D/V4z/k  (x  =  0.6328u,  r  =  1  km), 
for  various  values  o*  the  index  structure  constant 
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